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ABSTRACT 

This  paper  presents  a  negative  result  concerning  the  stable  evaluation  of  simplex  spline 
functions.  It  has  been  conjectured  that  a  great  deal  of  computational  effort  might  be  saved 
by  implementing  the  recurrence  relation  for  these  functions  in  a  clever  way.  The  main  result 
of  this  paper  is  that  there  is  no  clever  way  to  implement  the  recurrence  relation  once  the 
standard  recipe  for  constructing  spaces  of  simplex  spline  functions  has  been  followed. 
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SIGNIFICANCE  AND  EXPLANATION 


J  All  computational  schemes  for  simplex  spline  functions  to  date  rely  both  on  the  recur¬ 
rence  relation  for  these  functions  and  the  standard  construction  of  the  simplex  spline  basis. 
Under  these  conditions,  for  numerical  methods  for  computing  simplex  spline  functions  to 
be  as  useful  as  possible,  it  is  necessary  to  find  ways  of  implementing  the  recurrence  relation 
as  efficiently  as  possible.  This  paper  shows  that  “as  efficiently  as  possible”  is  still  not  very 
efficient.  This  implies  that  truly  fast  algorithms  (which  have  not  yet  been  developed)  will 
have  to  either  abandon  the  recurrence  relation  or  the  standard  construction  of  the  simplex 
spline  basis.  ~  •’  /?;  i' 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary  lies 
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THE  COMPUTATIONAL  COST  OF  SIMPLEX  SPLINE  FUNCTIONS 


Thomas  A.  Grandine1,2 

Since  de  Boor  defined  the  multivariate  simplex  spline  in  B76  .  many  of  its  mathe¬ 
matical  properties  have  been  studied  and  reported  (  B82  .  iD79],  M80  .  ...).  However, 
numerical  methods  for  evaluating  these  functions  have  not  been  studied.  The  goal  of  this 
paper  is  to  produce  an  algorithm  which  evaluates  simplex  spline  functions,  i.e.  linear 
combinations  of  simplex  splines.  This  can  be  done  in  a  straightforward  manner  by  using 
the  evaluation  algorithm  for  individual  simplex  splines,  described  in  G84  .,  to  evaluate 
all  non-zero  simplex  splines  at  any  point  x,  then  summing  up  the  values  using  the  spline 
function  coefficients  as  weights.  However,  since  this  algorithm  is  based  on  the  recurrence 
relation  for  simplex  splines,  M80),  it  seems  a  wasteful  approach.  In  particular,  many  of 
the  lower  degree  simplex  splines  which  need  to  be  evaluated  during  this  process  may  well 
be  evaluated  multiple  times.  It  seems  reasonable  (and  has  been  conjectured  by  Dahmen 
and  Micchelli  [DM81]  and  perhaps  others)  that  a  great  deal  of  computational  effort  might 
be  saved  by  taking  advantage  of  this  multiplicity.  This  paper  makes  an  attempt  to  do  just 
that. 

In  order  to  begin  this  attempt,  the  recurrence  relation  is  needed.  This  is  given  by 
M80 

n 

M{x\On.ol,...,on)  =  ^  a, M (*]<;,;„ P] . pt_,.pt  +  1, (1) 

1=0 
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Here  Pz  means  the  first  m  components  of  the  vector  z. 


The  second  necessary  step  is  to  clarify  just  which  simplex  splines  are  the  correct  ones 
to  deal  with  when  applying  this  recurrence  relation.  This  requires  the  construction  of  an 
appropriate  simplex  spline  basis.  Hollig  H82  and  Dahmen  and  Micchelli  DM82  have 
simultaneous!)  and  independent!)  given  versions  of  one  such  construction.  The  Hollig 
version  of  the  construction  is  summarized  here: 

To  construct  a  simplex  spline  basis  over  fi  j  Rn'.  start  with  a  triangulation  of  fi. 
A  triangulation  of  a  set  D  is  a  collection  T  of  sets  r  such  that 

(*)  voljr]  >  0 

(u)  .  L'rtf T  k]  =  Q 
(in')  |t|  fi  |r'j  =  \t  n  r'j,  t,t'  e  T 

(tv)  #r  =  m  +  1,  r  €  T, 

where  \A)  is  the  convex  hull  of  A  and  41 A  is  the  cardinality  of  A.  Thus,  each  [rj  is  a 
simplex  in  U.  Consider  the  slab  fl  •  :=  {(u,v)ju  €  H,r  t  Sfc},  where  E k  is  the 

standard  unit  simplex  in  RA  The  construction  boils  down  to  triangulating  this  slab  and 
using  the  vertices  of  each  of  the  resulting  simplices  as  knot  sets  for  the  simplex  splines. 
The  slab  can  be  triangulated  as  follows:  For  each  r  6  T,  consider  the  following  grid: 


Here,  tl}  is  intended  to  represent  ( tt ,  e} ),  where  t,  is  a  point  in  r,  and  e3  is  the  j-th  standard 
unit  vector  (with  eo  being  0).  Hence  tX]  is  a  vertex  of  rxEj.  (This  is  not  quite  true.  In  the 
actual  construction,  the  knots  are  pulled  apart.  This  means  that  ttJ  is  not  really  ( tt,ej ), 
but  rather  some  point  which  is  close  to  (/ , .  ) .  The  actual  details  are  quite  interesting  and 

useful,  but  not  of  importance  here).  Consider  all  non-decreasing  paths  from  tn,n  to  tm 
Such  a  path  is  denoted  by  a,,,  ...,  o„.  where  o,,  =  f,,,,,  on  —  ,  and.  if  Of  ~  ttJ,  then 

+  i  is  either  +  or  t +  The  set  of  all  paths  through  the  grid  for  t  for  all  t  £  T  give 
rise  to  a  simplex  spline  basis. 

Thus,  each  of  the  knots,  a,,  can  be  classified  according  to  where  it  was  found  along 
the  path  from  tn.o  to  fm,k-  Specifically,  given  ot  -  tl<3,  then  is  called  a  horizontal 
knot  if  and  Or j_i  =  t,  +  ]j.  If  Ot_\  =  ttj-i  and  ot+  j  =  tt)y+ j,  then  oi  is 

called  a  vertical  knot.  All  other  knots  are  called  corner  knots.  Note  that  these  names 
are  specific  to  a  given  simplex  spline,  i.e.  a  knot  which  is  of  one  type  in  one  simplex  spline 
may  be  of  a  different  type  in  a  different  simplex  spline.  This  classification  of  knots  has  no 
significance  except  to  relate  the  simplex  splines  to  the  paths  from  which  they  were  derived. 
The  paths  provide  the  key  to  piecing  lower  degree  simplex  splines  together,  as  will  be  made 
clear  shortly. 

Definition:  The  simplex  spline  Mt  M(-|oy.,...,ct,_i,oi+],  ...,o„)  is  said  to  con¬ 
tribute  to  the  simplex  spline  if,  in  (1),  A/,  necessarily  has  a  non-zero 

coefficient  corresponding  to  it. 

In  order  to  prove  the  upcoming  theorem  in  this  paper,  two  more  concepts  are  needed. 
The  first  notion  is  that  of  a  cut  corner.  A  path  with  a  cut  corner  is  obtained  from  a  path 
<To, ..., by  removing  a  corner  point.  This  amounts  to  finding  a  path  through  the  grid 


consisting  of  only  horizontal  and  vertical  moves,  except  for  one  place,  arbitrarily  chosen, 
where  a  diagonal  move  is  allowed.  The  second  notion  has  to  do  with  the  grids  themselves. 
In  addition  to  the  grid  for  r.  a  grid  may  also  be  considered  for  t  n  t'.  If  t  i  t' .  then  this 
grid  will  have  fewer  than  rn  -•  1  columns,  and  a  path  through  the  grid  will  only  produce  a 
simplex  spline  of  degree  k  -  c  -  m  1.  where  c  is  the  number  of  columns  of  the  grid. 

Theorem  1:  Assume  that  t,0  ?  t unless  i  i'  and  j  -  j'  ■  Then  the  only  simplex 
splines  of  degree  k  1  which  can  contribute  to  two  simplex  splines  of  degree  k  are  those 
which  conic  from  paths  with  cut  corners  and  those  which  come  from  paths  corresponding 
to  grids  for  t  r  r'.  where  *  ~  '  r'  -  m.  There  are  no  simplex  splines  of  degree  k  —  1  which 
contribute  to  three  or  more  simplex  splines  of  degree  k. 

Proof:  Take  any  simplex  spline  of  degree  fc,  say  A/,  deriving  from  a  path  correspond¬ 
ing  to  the  grid  for  t.  Then  the  only  simplex  splines  which  can  contribute  to  M  are  those 
obtained  by  deleting  a  knot,  say  along  that  path.  Suppose  o3  is  a  vertical  knot.  Then 
Mj  Af(-|<70, i,...,<t„)  can  only  contribute  to  Af,  since  the  only  way  to  get 
from  to  Oj4-i  along  a  path  is  through  Oj.  Suppose  a}  is  a  corner  knot.  Assume 

that  o}  -\  =  tj-c-i,*’.  Then,  since  is  a  corner  knot,  aJ  +  \  =  t3 Thus,  the  only 
way  to  get  from  i  to  c,.,  along  a  path  is  through  either  t;_£_i,£+i  or  one  of 

which  must  be  o}.  Thus.  A /,  can  contribute  to  at  most  two  simplex  splines  of  degree  d. 
Lastly,  suppose  o }  is  a  horizontal  knot.  Then  comes  from  a  path  along  a  grid  where 
the  column  corresponding  to  a ,  is  deleted.  Suppose  there  exists  t'  such  that  rOr' provides 
just  such  a  grid.  Then  contributes  both  to  A/  and  to  the  corresponding  simplex  spline 
from  the  grid  for  r\  but  to  no  others.  Otherwise,  M}  just  contributes  to  A/.  This  proves 


the  theorem. 


wv  v> ^  ^  V* *;■  -’.- v*  jtv  *■■'  r'wjv.  .-  Ji-.  «r;  ^;t 


/  V-T  .'  '7  ?  ?  ’  'J  -7  V 


A  computer  code  has  been  written  using  this  idea.  The  same  approach  to  implement¬ 
ing  the  recurrence  relation  that  is  employed  in  G84  is  used  here.  Here,  however,  it  is 
best  to  have  non-zero  as  few  of  the  a,  corresponding  to  vertical  knots  as  possible.  Unfor¬ 
tunately.  including  this  condition  results  in  a  problem  whose  computational  complexity  is 
equivalent  to  that  of  the  travelling  salesman  problem.  Since  the  goal  here  is  to  write  a 
last  code,  the  simple  heuristic  device  of  steering  the  method  toward  corner  and  horizontal 
knots  whenever  possible  is  employed.  This  does  not  behave  perfectly,  but  it  does  func¬ 
tion  adequately,  and  the  non-zero  o,  obtained  usually  correspond  to  corner  and  horizontal 
knots  when  such  solutions  are  possible.  Additionally,  the  code  keeps  track  of  which  lower 
degree  simplex  splines  have  already  been  computed  so  that  their  results  can  be  reused  if 
possible.  A  hashing  scheme  has  been  implemented  to  provide  nearly  direct  access  to  the 
values  of  these  already  computed  simplex  splines. 

Unfortunately,  the  code  behaves  quite  poorly  in  practice.  In  actual  comparison  with 
the  more  naive  approach,  i.e.,  the  one  in  which  no  intermediate  results  are  saved  and  all 
simplex  splines  are  computed  one  at  a  time,  the  naive  approach  is  the  clear  and  decisive 
winner.  The  clever  approach  takes  approximately  three  times  as  much  CPU  time,  and  this 
ratio  seems  not  to  change  with  either  the  degree  of  the  simplex  splines  or  the  complexity 
of  the  knots,  at  least  through  seventh  degree  polynomials  in  the  bivariate  case! 

The  explanation  for  this  lies  in  the  excessive  overhead  involved  in  keeping  track  of  the 
lower  degree  intermediate  results.  The  only  way  to  distinguish  between  values  of  different 
simplex  splines  is  to  somehow  distinguish  between  the  simplex  splines  themselves.  This 
can  only  be  done  by  enumerating  the  knots  which  make  up  the  simplex  spline.  Thus,  the 
amount  of  overhead  increases  with  the  degree  of  the  simplex  splines.  This  means  that  the 
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payoff  that  might  be  expected  by  making  the  degree  of  the  simplex  splines  sufficiently  large 
does  not  actually  occur. 

What's  worse.  Theorem  1  says,  in  some  sense,  that  this  cannot  be  improved  on.  While 
other  programmers  might  be  able  to  improve  the  code  sufficiently  to  make  it  competitive, 
it  cannot  be  done  in  such  a  way  that  the  performance  of  the  naive  method  can  be  bettered 
by  more  than  a  factor  of  2k~  ’.  Indeed,  such  a  program  would  have  to  have  some  overhead 
not  present  in  the  naive  code,  and  th  re  will  usually  be  lower  degree  simplex  splines 
which  contribute  to  only  one  higher  degree  simplex  spline  present.  Thus,  a  factor  of  2k~  1 
improvement  is  really  far  too  much  to  hope  for. 
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